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ABSTRACT 

Context. We present VLT/X-shooter spectroscopic observations in the wavelength rang c AA3000 - 23000A of the ex- 
tremely metal-deficient blue compact dwarf (BCD) galaxy PHL 293B containing a luminous blue variable (LBV) star 
and compare them with previous data. 

Aims. This BCD is one of the two lowest-metallicity galaxies where LBV stars were detected, allowing us to study the 
LBV phenomenon in the extremely low metallicity regime. 

Methods. We determine abundances of nitrogen, oxygen, neon, sulfur, argon, and iron by analyzing the fluxes of narrow 
components of the emission lines using empirical methods and study the properties of the LBV from the fluxes and 
widths of broad emission lines. 

Results. We derive an interstellar oxygen abundance of 12+log O/H = 7.71 ± 0.02, which is in agreement with pre- 
vious determinations. The observed fluxes of narrow Balmer, Paschen and Brackett hydrogen lines correspond to the 
theoretical recombination values after correction for extinction with a single value C(RP) = 0.225. This implies that 
the star-forming region observed in the optical range is the only source of ionisation and there is no additional source 
of ionisation that is seen in the NIR range but is hidden in the optical range. We detect three v = 1-0 vibrational lines 
of molecular hydrogen. Their flux ratios and non-detection of v = 2-1 and 3-1 emission lines suggest that collisional 
excitation is the main source producing H2 lines. For the LBV star in PHL 293B we find broad emission with P Cygni 
profiles in several Balmer hydrogen emission lines and for the first time in several Paschen hydrogen lines and in several 
He 1 emission lines, implying temporal evolution of the LBV on a time scale of 8 years. The Ha luminosity of the 
LBV star is by one order of magnitude higher than the one obtained for the LBV star in NGC 2363=Mrk 71 which 
has a slightly higher metallicity 12+logO/H = 7.87. The terminal velocity of the stellar wind in the low-metallicity 
LBV of PHL293B is high, ~ 800 km s~\ and is comparable to that seen in spectra of some extragalactic LBVs during 
outbursts. We find that the averaged terminal velocities derived from the Paschen and He I emission lines are by some 
~ 40 - 60 km s" 1 lower than those derived from the Balmer emission lines. This probably indicates the presence of the 
wind accelerating outward. 

Key words, galaxies: fundamental parameters - galaxies: starburst - galaxies: ISM - galaxies: abundances - stars: 
activity 



1. Introduction cellent laboratories to test metallicity effects on mass loss 

d ; (jlzotov et al. 1 120071) . 

Critical evolutionary phases in the life of massive st ars such 

as th e stage of luminous blue variables (LBVs; iConti I . The prevailing belief is that the efficiencies of stellar 

[1981 are hard to study due to the exceedingly short time winds in massive stars are significantly reduced at low 

scales involved. To enhance the chance of catching the LBV metallicities. However, there is growing evidence suggest- 

phenomenon. it is best to observe galaxies known to contain in § that massive stars do have stellar winds, even in the 

many massive stars. Blue compact dwarf (BCD) galaxies most metal-deficient BCDs known with interstellar oxygen 

are suitable objects for such studies as their star-forming abundances 12 -flog. O/H <, 7.6. j lzotov et al. | (jlMfl) and 

regions can harbor up to thousands of massive stars (e.g., iLegrand et al. | UMZD discovered a Wolf-Rayet (WR) stel- 

iThuan fc Izotovl [20051) . As the mass loss rate depends on lar population in I Zw 18 with an oxygen abundance 12+log 

the metallicity, BCDs of varying metal deficiency are ex- O/H = 7.17 ± 0.01 jlzotov et al. |, H999J). Furthermore, 

iThuan fc Izotov I (|1997l ) have observed stellar P Cygni pro- 
files for the Si iv A1394, A1403 absorption lines in the spec- 
Send offprint requests to: Y.I. Izotov, izotov@mao.kiev.ua trum f the emission-line galaxy SBS 0335-052E with an 
* Based on observations collected at the European Southern inter stellar oxygen abund ance of 12 + log O/H = 7.30 ± 
Observatory, Chile, ESO program 60.A-9442(A). q 01 (|I zo toy et al. I . [25571 ) . 
** The reduced data in Figures 1 and 2 are available at the 

CD S via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) o r Although we do not know the metallicities of these 

via |http: / /cdsweb.u-strasbg.fr / viz-bin/qcat? J / A+ A/ 533 / A25 stars, it appears that they were formed in the interstellar 
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medium with an extremely low metallicity and do possess 
stellar winds. The properties of these stellar winds may sig- 
nificantly diffe r from those of their high- metallicity coun- 
terpa rts (e.g., iGuseva et all 120001: ICrowther fc Hadfieldl . 
l2006h . 

To find more obje cts with obse r vation al signatures of 
massive star outflows, llzotov et al. I (|2007l ) have assembled 
a sample of about 40 emission-line dwarf galaxies which 
exhibit broad components in their strong emission lines, 
mainly in H/3, [O in] AA 4959, 5007, and Ha. Except for 
four objects which appe ar to contain massive black holes 
(jlzotov fc Thuan I . I2008I ). the broad emission of all the other 
objects can be attributed to some evolutionary stage of 
massive stars and to their interaction with the circumstellar 
and interstellar medium, e.g. WR star s, supernovae, super - 
bubbles or LBV stars. The attention of llzotov et al. I (|2007l) 
was drawn to the spectrum of the BCD PHL 293B = SDSS 
J2230-0006, with 12 + log O/H = 7.72, which shows broad 
B aimer hydrogen emission lines with P Cygni profiles, that 
are char acteristic of a LBV sta r. Several months after this 
finding, iPustilnik et al. I (|2008l ) discovered another bright 
LBV in the BCD DDO 68 with interstellar 12 + log O/H 
= 7.14. The study of thes e newly discovered LBV stars was 
the focus of the paper by llzotov fc Thuanl (|2009h . 

The LBV phase represents a transition in the late evo- 
lution of all stars with initial masses greater than about 
50 Mq from the stage of O stars burning hydrogen on 
the main-sequence to that of core-helium burning WR 
stars. During this phase, the stars lose large amounts of 
mass in recurrent explosive events. In the HR diagram, 
LBVs lie just to the left o f the H umphreys-Davidson limit 
(|Humphrevs fc Davidson I 1 1994T ) . beyond which no stars 
are observed. They define the locus of an instability that 
prevents fur ther redward evolution, e.g. the "strange mode 
instability" (jKiriakidis et al. l . fl993T) . 

The exact mechani sm giving rise to a LBV outburst 
is sti ll unknown (e.g., iMaeder fc MevneFI . [2000; Ma ssev I . 
1200311 . even though many LBVs at sol ar metallicity have 

. The situation is 
(|2011| ) noted that 
arc otten detected in spiral galaxies but rarely 
in dwarf galaxies and the metallicities of most of these 
dwarf galaxies are still unknown. Exceptions are the Small 
Magellanic Cloud with 12+log O/H - 8.0 and the dwarf 
cometary galaxy NGC 2366 with 12+log O/H = 7.87. In the 
latter galaxy, the LBV star has been investiga ted in detail 
(jDrissen et al. I . I1997L 120011: iPetit et al. 1 . 120061) . Repeating 
high-quality spectroscopic observations of LBV stars in 
more extremely metal-deficient galaxies should help to bet- 
ter understand how the LBV phase depends on metallicity. 

We present here new VLT/X-shooter spectroscopic ob- 
servations of PHL 293B over a larger wavelength range of ~ 
A3000-23000A as compared to the SDSS (~ A3800-9200A) 
an d VLT/UVES (~ A3000-7000A) optical data discussed 
bv llzotov fc Thuanl (|2009f ). These' new observations allow 
us to derive element abundances in the H II region, to study 
whether any hidden star formation is present in this galaxy, 
and to obtain emission line parameters of the LBV star for 
the epoch ~ 8 years after previous observations. The new 
observations are described in Sect. [5J In Sect. [3] we discuss 
the properties of the host galaxy. We derive the element 
abundances, discuss extinction and hidden star formation 
and excitation mechanisms of the H2 emission in the NIR 
range. The properties of the LBV broad emission are anal- 
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2. Observations and data reduction 

A new spectrum of PHL 293B was obtained using the 8.2 
m Very Large Telescope (VLT) on 2009 August 16 [ESO 
program 60.A-9442(A)]. The observations were performed 
at an airmass ~ 1.4 using the three-arm echelle X-shooter 
spectrograph mounted at the UT2 Cassegrain focus. In the 
UVB arm with wavelength range A3000 - 5595A, four equal 
exposures of 900s with a l"xll" slit were obtained result- 
ing in a total exposure time of 3600s. The angular scale of 
1" at a distance of 22.7 Mpc to PHL 293B taken from the 
NASA/IPAC Extragalactic Database (NED)Q corresponds 
to a linear scale of 110 pc. The distance of PHL 293B 
has been obtained from the radial velocity corrected for 
Virgo Infall with a Hubble constant of 73 km s -1 Mpc -1 . 
Nodding along the slit was performed according to the 
scheme ABBA with the object positions A or B differ- 
ing by 5" along the slit. In the VIS arm with wavelength 
range A5595 - 10240A, eight equal exposures of 450s with a 
0'.'9xll" slit were obtained requiring a total exposure time 
of 3600s. Nodding along the slit was performed according to 
the scheme AABBBBAA. In the NIR arm with wavelength 
range A10240 - 24800A, twelve equal exposures of 60s with 
a 0'.'9xll" slit were obtained resulting in a total exposure 
time of 720s. Nodding along the slit was performed accord- 
ing to the scheme AAABBBBBBAAA. The above instru- 
mental set-up resulted in resolving powers A/AA of 5100, 
8800 and 5100 for the UVB, VIS and NIR arms, respec- 
tively. The seeing was ~ 1". 

The Kitt Peak IRS spectrophotometric standard star 
Fcige 110 was observed on 2009 August 15 with a slit of 
5"xll" for flux calibration and several bright stars with 
the same slits as those used for PHL 293B were observed 
on 2009 August 15 to correct for telluric absorption. We 
use the spectrum of Feigc 110 for the flux calibration of the 
spectra obtained with the UVB and VIS arms. Since this 
star was observed with a wide slit, no light losses need to 
be taken into account for the flux calibration. As for the 
NIR arm, we use the spectrum of one of the telluric stars, 
Hip 15378, of spectral type A0V with a V magnitude of 
7.86 mag, to perform the flux calibration. Since there is 
no available absolute spectral distribution for that star, wc 
adopted the spectral energy distribution of another A0V 
star, Vega, with a V magnitude of 0.0 mag and scaled it to 
the V magnitude of Hip 15378. Spectra of thorium-argon 
(Th-Ar) comparison arcs were obtained and used for wave- 
length calibration of the UVB and VIS arm observations. 
For the wavelength calibration of the NIR spectrum we use 
night sky emission lines. 

The two-dimensional UVB and VIS spectra were bias 
subtracted and flat-field corrected using IRAF0. The two- 
dimensional NIR spectra were corrected for a dark current 



1 NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the National Aeronautics and 
Space Administration. 

2 IRAF is the Image Reduction and Analysis Facility dis- 
tributed by the National Optical Astronomy Observatory, which 
is operated by the Association of Universities for Research 
in Astronomy (AURA) under cooperative agreement with the 
National Science Foundation (NSF). 



Y.I.Izotov et al: VLT/X-shooter observations of the low-metallicity BCD PHL 293B 



3 



10 



I 



Ml 





3000 



PHL 293B 




3200 



3400 



3600 



3800 



4000 



4200 



4400 




10 



10 



10 



4600 



4800 



5000 



5200 



5400 



6000 



6200 



6400 



6600 



6800 



5600 



5800 



7000 



7200 



8800 9000 9200 9400 9600 

Wavelength (A) 



9800 



10* 



1 1 1 1 

o o 

i-i CQ CO 

m n n 

OJ i> o 
IN 

: §la 


- [Arlll] 7751 

' Pa21 B374 

■ Pa20 B392 
Pal9 8413 

■ PalB 8438 
' Pal7 B470 

• Pal6 8505 
' Pal5 B54B 

• Pal4 8601 - 

- Pal3 B6B7 



1,1,1,1,1,1,1, 


7400 7600 


7800 8000 8200 8400 8600 


1 I 1 I 

o 

. Is Co in N £ 
co in co ^2 

CD I CO O 02 
■ CO CO CO OO , , 

" co w -h ok 
- ^ ^ ^ ^ co 
.co co co co — 1 

. Ph Ph Ph Ph 

" L . _ jl_ 1 


1 I 1 I 1 

CO 
CO 

~h CO 

co <s 

OJ , — 

oo fji 
CO i — 

Ph 


1 ' 1 ' 1 ' 

cj m 

o 

in o 

ID 'O ~ 

cd to 

Ph Ph 

i. JL : 


1,1,1,1,1,1,1, 



Fig. 1. Flux-calibrated and rcdshift-corrected VLT/X-shooter UV and optical spectrum of PHL 293B. Note the broad 
emission with P Cygni profiles in the hydrogen lines US A4101 through Ha A6563. 



and divided by the flat frames to correct for the pixel sensi- 
tivity variations. Cosmic ray hits of all UVB, VIS and NIR 
spectra were removed using routine CRMEDIAN. The re- 
maining hits were later removed manually after background 
subtraction. For each of the UVB, VIS and NIR arms we 
separately coadded spectra with the object at the position 
A and spectra with the object at the position B. Then, 
the coadded spectrum at the position B was subtracted 
from the coadded spectrum at the position A. This re- 
sulted in a frame with subtracted background. We used 
the IRAF software routines IDENTIFY, REIDENTIFY, 
FITCOORD, and TRANSFORM to perform wavelength 
calibration and correct for distortion and tilt of each 
frame. The one-dimensional wavelength-calibrated spectra 
were then extracted from the two-dimensional frames us- 
ing the APALL routine. We adopted extraction apertures 
of l"x2'/5, 0'.'9x2'.'5 and 0'.'9x2'.'5 for the UVB, VIS and 
NIR spectra, respectively. 

The resulting flux-calibrated and rcdshift-corrected 
UVB and VIS spectra of PHL 293B are shown in Fig. [1] 
and the resulting flux-calibrated and redshift-corrected NIR 
spectrum is shown in Fig. [2j Despite the fact that UVB and 
VIS spectra, on the one hand, and the NIR spectrum, on 
the other hand, were flux-calibrated with the use of two 
different stars observed with different apertures, the levels 
of the continuum at ~ 10000A in the VIS and NIR spectra 
agree to better than 10%. Therefore, no adjustment of the 
VIS and NIR spectra was applied. Strong broad hydrogen 
emission lines with P Cygni profiles are present in the UVB 
and VIS spectra, similar to the ones seen in the UVES spec- 



trum (jlzotov fc Thuanl . l2009f) . Two hydrogen lines in the 
NIR spectrum (Fig. Pa/3 and Pa7, also show the broad 
components. 

3. The properties of the host galaxy 

3.1. Element abundances 

We derived element abundances from the narrow emission- 
line fluxes, using a classical semi-empirical method. These 
lines trace the interstellar medium (ISM) of PHL 293B. 
The fluxes in all spectra were measured using Gaussian 
fitting with the IRAF SPLOT routine. The lcr errors of 
the line fluxes were calculated from the data in the non- 
flux-calibrated spectra using SPLOT. The method is based 
on Monte-Carlo simulations, in which random gaussian 
noise with the dispersion obtained from the continuum 
near the line is added to the noise-free spectrum. To ob- 
tain good error estimates we have run 200 simulations per 
emission line. These errors were propagated in the calcu- 
lations of the elemental abundance errors. The fluxes were 
co rrected for bot h extin ction, using the reddening curve 
of ICardelli et al. I (|1989t) . and underlying hydrogen stellar 
absorption, derive d simultaneou s ly by an iterative proce- 
dure described by llzotov et al. I (119941 ) and using the ob- 
served decrements of the narrow hydrogen Balmer lines 
in the UVB and VIS spectra. The extinction coefficient 
C(H/3) and equivalent width of hydrogen absorption lines 
EW(abs) are derived in such a way as to obtain the clos- 
est agreement between the extinction-corrected and theo- 
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Wavelength (yum) 

Fig. 2. The same as in Fig. Q] but for the near-infrared (NIR) range. Note the broad emission in the Pa/3 and Pa7 lines. 
Broad emission in the Paa line is also present but it is strongly absorbed by telluric lines. 



retical recombination Balmcr hydrogen emission-line fluxes 
normalised to the H/3 flux. It is assumed that EW(abs) is 
the same for all hydrogen lines. This assumption is justi- 
fie d by the evolutionary ste l lar po pulation synthesis models 
of iGonzalez Delgado et all (|2005ft . The derived C(H/3) and 
EW a b s are applied for correction of the emission-line fluxes 
in the entire wavelength range AA 3000 - 22650A. 

The extinction-corrected total fluxes 100 x /(A) //(H/3) 
of the narrow lines, and the extinction coefficient C(H/3), 
the equivalent width of the H/3 emission line EW(H/3), the 
H/3 observed flux F(H/3), and the equivalent width of the 
underlying hydrogen absorption lines EW(abs) are given 
in Table Q] The extinction-corrected flu xes are generally 
in goo d agreement with those derived by llzotov fc Thuanl 
(12009ft . 

The physical conditions, and the ionic and total heavy 
clement abundanc es in the H I I regio n of PHL 293B were 
derived following llzotov et all (|2006h (Tabic [2]). In par- 
ticular for 2+ , Ne 2+ , and Ar 3+ , we adopt the tem- 
perature T^ ir (0 in) directly derived from the [O in] 
A4363/(A4959 + A5007) emission-line ratio. The electron 
temperatures T e (0 n) an d T P (S ill) w e re de rived from 
the empirical relations by llzotov et al.l (|2006l) based on 
photoionisation models of H II regions. By direct meth- 
ods we derived also Tf r (0 n)= 15820±1460 K from the 
[O n]A(3726+3729)/A(7320+7330) emission-line ratio and 
Tf r (S lll)= 15640±3230 K from the [S m]A6312/A9069 
emission-line ratio. Despite their higher errors these tem- 
peratures are in agreement with the temperatures T e (0 il) 
and T e (S in) derived from the empirical relations. We used 



T e (0 n) for the calculation of 0+, N+, S+, and Fe 2+ 
abundances and T e (S III) for the calculation of S 2+ and 
Ar 2+ abundances. The electron number densities N e (0 II), 
N e (S II) and N e (Ar iv) were obtained from the [O n] 
A3726/A3729, [S n] A6717/A6731 and [Ar iv] A4711/A4740 
emission-line ratios, respectively. For these low electron 
number densities (Table [2]), the collisional deactivation of 
the upper levels of forbidden optical and NIR emission lines 
is negligible as compared to the spontaneous transitions. 
The element abundances then do not depend sensitively 
on N e and the low-density limit holds for the abundance 
determination in PHL 293B. The ionic and total O, N, 
Ne, S, Ar, and Fe abundances derived from the forbidden 
emission lines are given in Table [5] The oxygen abundance 
12+log O/H = 7.71±0.02 and abundance ratios log N/O = 
-1.59±0.05 and log Ne/O = -0.77±0.03 are consistent with 
the respecti ve values 7.72±0 1, -1 55±0.02 and -0.74±0.02 
derived by llzotov fc Thuanl (|2009t ). All abundance ratios 
in Table [2] are typical for th e low-metallicity emission-line 
galaxies (jlzotov et al.l . l2006h . 



3.2. Extinction and hidden star formation 

In al l previous studies o f BCDs except for those of Mrk 
59 in llzotov et all (12009ft . II Zw 40 Mrk 71, Mrk 996 and 
SBS 0335-052E in llzotov fc Thuanl (120111V the NIR spec- 
tra of BCDs were obtained in separate JHK observations, 
and there was no wavelength overlap between the optical 
and NIR spectra. The fact that the spectrum of PHL 293B 
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Table 1. Narrow emission- line fluxes with la errors 



Line 


100xF(A)/F(H,fl) a 


100x/(A)//(H/3) b 


Line 




100xF(A)/F(H/3) a 


100x/(A)//(H^) ^, 


3f88 He I 


2.98±0.22 


3.73±0.35 


6717 


"TTt T 

[S ii] 


5.28L0.08 


4.45L0.28 


3683 H20 


0.55±0.04 


2.39±0.28 


6731 


[S 11] 


3.96±0.07 


3.33±0.24 


3687 Hf9 


0.47±0.04 


2.29±0.30 


7065 


He 1 


2.66L0.05 


2.18±0.19 


3692 Hf8 


1.02±0.04 


2.92±0.28 


7136 


[Ar III] 


4.75±0.07 


3.89L0.26 
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0.98±0.04 


2.90±0.27 


7281 


He 1 


0.66L0.02 


0.54±0.10 


3704 Hf6 
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2.90±0.28 


7320 


[O 11] 


1.07L0.03 


0.87±0.12 


3705 He I 


0.52±0.05 


0.61±0.12 


7330 
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0.84±0.12 


3712 H15 


1.21±0.04 


2.90±0.27 


7751 


[Ar 111] 


1.42±0.03 
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3722 H14 
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4.56±0.34 


8392 


Pa20 


0.35±0.02 
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3726 [0 n] 


18.42±0.27 


21.43±0.74 


8413 


Pal9 


0.32L0.02 


0.53±0.09 


3729 [O n] 


25.56±0.38 
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8438 


Pal8 
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3734 H13 
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8446 


O 1 
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3750 H12 


2.18±0.05 


4.23±0.32 


8470 


Pal7 
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3770 Hll 
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8505 
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0.76L0.11 


3797 H10 


4.08±0.09 


6.39±0.40 


8548 


Pal5 


0.65L0.02 


0.77L0.11 


3820 He I 
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1.02±0.16 


8601 


Pal4 


0.90L0.03 
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3835 H9 


5.91±0.11 


8.41±0.45 
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Pal3 


0.97±0.02 


0.99±0.13 


3868 [Ne in] 


37.69±0.55 


43.18±1.08 


8753 


Pal2 


1.28±0.04 


1.23±0.14 


3889 He I+H8 


17.35±0.26 


21.59±0.74 
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Pall 


1.71L0.04 


1.54±0.16 
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0.55±0.05 
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9017 
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3970 H7 
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4068 [S n] 
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0.18±0.07 
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4340 H7 
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47.39±1.09 
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27.38±0.49 
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0.83L0.09 


0.60±0.10 
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0.45L0.04 


0.46±0.10 


17006 He 1 


0.40±0.09 


0.26±0.07 


4740 [Ar iv] 


1.32±0.04 


1.33±0.17 


17367 BrlO 


1.23±0.08 


0.83L0.11 


4ool rip 


1 nn nn_Li a o 


1 AA AA_1_1 A A 


18179 Br9 c 


0.80L0.08 


0.51L0.09 


4921 He 1 


1 02+0 03 

1 ■ \ > — — 1 — \ 1 • \ f r J 


1.00±0.14 


18756 Paa c 


2.57±0.10 


1.65±0.15 


4959 [O in] 


188.62±2.68 


184.30±2.82 


19451 Br<5 c 


0.64±0.10 


0.43±0.09 


4986 [Fe 111] 


0.33±0.02 


0.32±0.08 


20337 H 2 1-0 S(2) 


0.27L0.05 


0.17L0.05 


5007 [O in] 


558.58±7.93 


542.49±7.85 


20586 He i c 


0.66±0.06 


0.42±0.08 


5016 He 1 


1.88±0.04 


1.82±0.19 


21220 H 2 1-0 S(l) 


0.67±0.06 


0.42±0.08 


5876 He 1 


9.79±0.15 


8.74±0.41 


21661 Br 7 


3.62L0.11 


2.28±0.18 


6300 [O 1] 


1.14±0.03 


0.99±0.13 


22195 [Fe ill] 


0.35±0.09 


0.22±0.07 


6312 [S in] 


1.47±0.04 


1.27±0.15 


22230 H 2 1-0 S(0) 


0.34L0.10 


0.21±0.07 


6363 [O 1] 


0.62±0.07 


0.53±0.10 


c{W) 


0.225 


6563 Ha 


333.33±4.73 


284.26±4.45 


EW(H/3) d 


112 




6583 [N 11] 


2.04±0.05 


1.73±0.18 


F(R/3) e 


83.8 




6678 He 1 


3.02±0.06 


2.56±0.21 


EW(abs) d 


1.15 





a observed flux. 

b flux corrected for extinction and underlying stellar absorption. 
c affected by telluric absorption. 

d equivalent width of hydrogen absorption lines in A. 
e in units 10~ 16 erg s -1 cm -2 . 



has been obtained simultaneously over the entire optical 
and near-infrared wavelength ranges, avoids adjusting un- 
certainties and permits us to compare directly the optical 
and NIR extinctions. 

The agreement between the corrected values of the 
hydrogen line fluxes in both the optical and NIR 
ranges (Table U) a n d the oretical recombination values by 
iHummer fe Storevl (|l987f) implies that a single C(H/3) 
[or A(V)j can be used over the whole 3000 - 22650A 
range to correct line fluxes for extinction. In particu- 



lar, the extinction-corrected flux of the Br 7 emission 
line J(Br7)//(H/3) = 0.2 3±0.02 is in agreement with the 
IHummer fe Storevl (|1987l ) recombination value of 0.25 for 
the electron temperature T e = 15000K and the electron 
number density N e = 100 cm~ 3 . The fact that A(V) does 
not increase when going from the optical to the NIR wave- 
length ranges implies that the NIR emission lines do not 
probe more extinct regions with hidden star formation as 
compared to the optical emission lines. This appears to be 
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Table 2. Physical conditions and element abundances 



Table 3. H2 emission-line fluxes relative to the H 2 2.122/im 
flux 



Property 
Tf (O in), K 
T e (0 n), K 
Tf r (0 11), K 

T e (S III), K 

Tf r (S m), K 
iV e (0 11), cm- 3 
N e (S 11), cm" 3 
iV e (Ar iv), cm~ 3 

0+/H+, (xlO 5 ) 
2+ /H+, (xlO 5 ) 
3+ /H+, (xlO 6 ) 
O/H, (xlO 5 ) 
12+log O/H 

N+/H+, (xlO 6 ) 
ICF(N) a 
N/H, (xlO 6 ) 
log N/O 

Ne 2+ /H+, (xlO 5 ) 
ICF(Ne) a 
Ne/H, (xlO 5 ) 
log Ne/O 

S+/H+, (xlO 6 ) 
S 2+ /H+, (xlO 6 ) 
ICF(S) a 
S/H, (xlO 6 ) 
log S/O 

Ar 2+ /H+, (xlO 7 ) 
Ar 3+ /H+, (xlO 7 ) 
ICF(Ar) a 
Ar/H, (xlO 7 ) 
log Ar/O 

Fe 2+ /H+, ( xlO 6 ) (4658) 

ICF(Fe) a 

Fe/H, (xl0 6 )(4658) 

log Fe/O (4658) 

Fe 2+ /H+, ( xlO 6 ) (4986) 

ICF(Fe) a 

Fe/H, (xl0 6 )(4986) 

log Fe/O (4986) 



Value 



16690±360 
15220±300 
15820±1460 
15550±300 
15640±3230 
69±47 
76±166 
1680±530 

0.45±0.03 
4.63±0.25 
1.02±0.15 
5.18±0.25 
7.71±0.02 

0.12±0.01 

10.8 
1.34±0.12 
-1.59±0.05 

0.86±0.05 

1.04 
0.89±0.06 
-0.77±0.03 

0.07±0.01 
0.59±0.08 

2.30 
1.52±0.17 
-1.53±0.05 

1.47±0.10 
1.03±0.13 

2.22 
3.25±0.38 
-2.20±0.06 

0.05±0.01 

16.0 
0.76±0.22 
-1.83±0.13 

0.06±0.01 

16.0 
0.93±0.24 
-1.75±0.11 



Ionisation correction factor. 



a general result fo r BCDs (e.g. IVanzi et al. I . l2000l l2002t 
llzotov et all l2009tlkotov fe Thuanl . l201lh . 



3.3. H-2 emission 

Molecular hydrogen lines do not originate in the H II re- 
gion, but in neutral molecular clouds. In the near-infrared, 
they are excited through two mechanisms. The first one is a 
thermal mechanism consisting of collisions between neutral 
species (e.g., H, H 2 ) resulting from large-scale shocks due 
to the interaction of stars with molecular clouds or/and 
cloud-cloud collisions. The second one is the fluorescent 
mechanism due to the absorption of ultraviolet photons 
from the hot stars. By comparing the observed line ratios 



Line 




Observations 


Model" 








Fluor. 


Collis. 


1.314 H 2 


3-1 Q(l) 


<0.2 


0.6 


0.0 


2.034 H 2 


1-0 S(2) 


0.4±0.2 


0.5 


0.3 


2.122 H 2 


1-0 S(l) 


1.0 


1.0 


1.0 


2.223 H 2 


1-0 S(0) 


0.4±0.2 


0.6 


0.3 


2.248 H 2 


2-1 S(l) 


<0.3 


0.5 


0.0 



adopt their model 1 for fluorescent lines and model SI for colli- 
sionally excited lines. 



with those predicted by mode ls, such as those calculated by 
iBlack fc van Dishoeck I (|1987T) , it is possible to discriminate 
between the two processes. In particular, line emission from 
vibrational levels v>2 is virtually absent in collisionally ex- 
cited spectra, while they are relatively strong in fluorescent 
spectra. 

Three H 2 lines, 2.034 ^m 1-0 S(2), 2.122 fim 1-0 S(l) 
and 2.223 /im 1-0 S(0), are detected in the NIR spectrum of 
PHL 293B. These lines are labelled in Fig.[5]and their fluxes 
relative to that of the strongest 2.122 (im 1-0 S(l) line with 
the la errors are shown in Table [31 The two last columns of 
the Ta ble show theoretical ratios bv lBlack fc van Dishoeck 1 
(119871 ) for fluorescent and collisional excitation, respec- 
tively. It is seen from the Table that the observed line flux 
ratios of the lines from the vibrational level v = 1 are in 
agreement with the values expected for both the fluorescent 
and collisional excitation. However, no emission lines from 
the vibrational levels v > 2 were detected while they should 
be present if the fluorescent excitation is important. 

We show in Table |3j the la upper limits of the two 
strongest H 2 emission lines, 1.314 /jm 3-1 Q(l) and 2.248 
^m 2-1 S(l), originating from the upper levels v >2 and 
label their location in the spectrum on Fig. [2] In the case 
of excitation dominated by fluorescence these lines would 
be as strong as 2.034 /xm 1-0 S(2) and 2.223 jum 1-0 S(0) 
lines. Therefore, we conclude that H 2 vibrational states in 
the spectrum of PHL 293B arc mainly collisio nally excited. 
This i s at variance with co nclusions mad e by Vanzi et all 
(12000H for SBS 0335-05 2E, IVanzi et al. I (12008H for II Zw 
40. llzotov et al. I (|2009t) for Mrk 59, and llzotov fc Thuanl 



(|201lD for II Zw 40, Mrk 71, Mrk 930, Mrk 996 and SBS 
0335-052E. There it has been found that fluorescence is 
the dominating excitation mechanism of H 2 lines in BCDs. 
Apparently, the difference is related to the low H/3 luminos- 
ity of PHL293B that is ~ 6 times lower than that in Mrk 71 
and more than one orde r of magnitude lowe r than that in 
other B CDs studied bv IVanzi etall (|2000D. IVanzi et al I 
(|2008D . llzotov et all (|2009l ) and llzotov fc Thuanl 1)201 lh . 
indicating a low intensity of the UV radiation in PHL293B. 



4. Broad emission from the LBV star 

The detection of the LBV star in PHL 293B indicates 
that high-mass stars are present in its star forming re- 
gion. The equivalent width of the narrow H/3 emission line 
EW(H/3) = 112A in PHL 293B (Table HJ corresponds to a 
starburst age of ~ 4.5 M yr for the PHL 293B metallicity 
(|Schaerer fc Vaccal Il998l ). At this age ma ssive stars with 
mass es up to ~ 60 M still can be present dMevnet et al. I . 
1199411 and one of them was likely evolved to the LBV star 
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we see today. Below we consider the properties of the LBV 
broad emission lines seen in the VLT/X-shooter spectrum. 

4.1. Fluxes 

In Tabic 0] we show the observed fluxes of the broad and 
narrow components of hydrogen H7, H/3 and Ha emission 
lines for three epochs encompassing ~ 8 years. For the SDSS 
and UVES data, we adopte d fluxes and respective errors 
from llzotov fe Thuanl (|2009f ). For the X-shooter data, us- 
ing the IRAF/SPLOT routine, we decomposed the H7, H/3 
and Ha lines into a narrow emission gaussian component, 
a broad emission gaussian component and an absorption 
gaussian component. An example of the line decomposition 
is shown in Fig. [3] for the H/3 line profile. It is seen from 
the Figure that the fit nicely reproduces the observed line 
profile. Fluxes of the broad lines in Table |4] are those corre- 
sponding to the fluxes of broad gaussians with the respec- 
tive errors which are derived similarly to those for narrow 
lines according to prescriptions in Sect. 13.11 We note that 
real errors of the fluxes can be significantly larger due to, 
e.g., uncertainties of the flat- field correction, background 
subtraction, telluric absorption, flux calibration, which are 
difficult to be accounted for. 

We also note for the first time the presence of 
broad components of Pa7 and Pa/3 emission lines in 
the VLT/X-shooter spectrum with the observed fluxes of 
(4.7±0.7)xKT 16 and (8.1±0.5) x lfT 16 erg s" 1 cm" 2 , re- 
spectively. 

A comparison of fluxes in Table Q] for different epochs 
shows that they vary by Sa 50% and are the lowest for the 
VLT/UVES observations for both the broad and narrow 
components. The fluxes of broad components in SDSS and 
VLT/X-shooter spectra are in agreement, but fluxes of nar- 
row components are higher in the SDSS spectrum. The nar- 
row lines are most likely produced in the H 11 region ionised 
by many O stars. Therefore, we would not expect temporal 
variations of these lines. Probably, the observational un- 
certainties and different extraction apertures play a role. 
The observational uncertainties may be caused by varying 
weather conditions, uncertainties in acquisition, and dif- 
ferent epochs of observation of the object and the stan- 
dard star. The larger fluxes of narrow emission lines by 
~25% in the SDSS spectrum can be explained by a larger 
round 3" aperture as compared to the smaller extraction 
aperture for the VLT/X-shooter spectrum. The different 
apertures, however, should not influence the fluxes of the 
broad lines because these lines are emitted by the com- 
pact expanding envelope a round the LBV s t ar. A s for the 
VLT/UVES observations, llzotov fc Thuanl (|2009f) do not 
present extraction apertures. These observations were done 
under non-photomctric conditions, according to the ambi- 
ent conditions database of the Paranal observatory. This 
may be the reason why the fluxes of the emission lines in 
the VLT/UVES spectrum are the lowest. 

We also present in Table [4] the broad-to-narrow line flux 
ratios. At variance to the line fluxes these flux ratios should 
not be sensitive to the uncertainties in the flux calibra- 
tion, although they depend on the extraction apertures. 
The largest differences of the flux ratios are for the Ha 
line with the largest value in the X-shooter observations. 
If, however, the aperture correction by a factor of ~1.25 to 
the narrow Ha line flux in the X-shooter observations is ap- 
plied, the broad-to-narrow Ha line flux ratios in the SDSS 



1 

- 






- 


J 

spectrum and fit / _ " 


r 










s 

X 




gaussians ^ 
















\ / 








residua 1 

















4860 4880 4900 4920 

Wavelength (A) 



Fig. 3. Decomposition of the H/3 profile into absorption 
and narrow and broad emission gaussian profiles. The ob- 
served spectrum and the fit to it are shown by solid lines, 
gaussian profiles arc shown by dashed lines, and the residual 
spectrum is shown by the dotted line. For a better view, the 
spectrum, the fit and gaussians are shifted up by 8xI0~ 17 
erg s _1 cm _2 A" 1 . 



and X-shooter data would be in good agreement. The same 
aperture correction for the UVES data results in a Ha flux 
ratio, which is ~40% lower. Apparently, some variations of 
the broad emission-line fluxes in PHL 293B on a time scale 
of ~ 8 years were present. These temporal variations are 
followed by variations of the broad line widths, discussed 
below in Sect. H21 

The broad Ha line luminosity of the LBV in PHL 293B 
derived from the X-shooter spectrum and corrected for the 
extinction with C(H/3) = 0.225 is 9.0 x 10 38 erg s _1 . It com- 
pares well with the broad Ha luminosity of ~ f 39 erg s _1 
deri ved from the SDSS s pectrum adopting C(H/3) = 0.245 
from llzotov et al."1 (|2007t) . Note that the broad Ha luminos- 
ity of the LBV in PHL 293B of 6.0 x 10 38 erg s" 1 , reported 
bv llzotov fc Thuanl (|2009f ). was not corrected for extinction. 
Thus, the Ha luminosity of the LBV in PHL 293B is one 
order of magnitude higher than L(Ha) ~ 10 38 erg s _1 of 
the VI star in NGC 2363, th e third lowest-metallicity LBV 
known (jDrissen et al. I . l200lh after those in PHL 293B and 
DDO 68, and two orders of magnitude higher than L(Ha) 
- 9.4xI0 36 erg s" 1 of the LBV in DDO 68. We also can 
estimate the brightness of the LBV in the continuum. The 
intensity of, e.g., H7 and H/3 at the bottom of the blue 
absorption features in the VLT and SDSS spectra of PHL 
293B is 0.8 that of the continuum. Since the light of the 
LBV star is diluted by the emission of the galaxy and as- 
suming that the intensities at the bottom of these absorp- 
tion f eatures in the spectrum of the LBV itself are zero 
(e.g., iKudritzki fc PulsT l2000f) . we adopt that the fraction 
of the LBV continuum in the total continuum of PHL 293B 
is <~ 20%. With this fraction and using the apparent SDSS g 
magnitude of 18.2 inside the 3" fiber aperture of the SDSS 
spectrum we obtain an apparent magnitude for the LBV of 
g = 19.9, which is not corrected for extincion. This corre- 
sponds to an absolute magnitude M g = -f f .9 at a distance 
D = 22.7 Mpc. If the extinction coefficient C = 0.225 is 
taken into account, then the LBV star would be ~ 0.5 mag 
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Table 4. Observed fluxes of the broad and narrow hydrogen emission lines at different epochs 



Date 


Broad line flux" 


Narrow line flux™ 


Broad-to- 


narrow line fj 


ux ratio 




H 7 


Up 


Ha 


H 7 


Hp 


Ha 


H 7 


up 


Ha 


2001.08.22" 




26.6±4.5 


97.8±6.3 


46.1±2.1 


104.8±4.6 


351.0±7.0 




0.25±0.04 


0.28±0.02 


2002.11. 08 c 


7.5±0.2 


16.7±0.3 


55.2±0.6 


32.6±0.5 


71.6±0.8 


213.6±2.1 


0.23±0.01 


0.23±0.01 


0.26±0.01 


2009.08.16 d 


7.4±0.2 


22.8±0.3 


97.0±0.7 


36.0±0.5 


83.8±1.2 


279.3±4.0 


0.21±0.01 


0.27±0.01 


0.35±0.01 



a in units of 10~ 16 erg s _1 cm -2 . 

"fluxes are from the SDSS spectrum obtained with the round aperture of 3" in diameter (|Izotov fc T huan. 2003). 

c fluxes from the VLT/UVES spectrum obtained with a 1" wide slit under non-photometric conditions (|Izoto*v fc Thuanl. [20091 ). 
d fluxes from the VLT/Xshooter spectra within the apertures l"x2 / /5 for the H 7 and H/3 lines and 0'.'9x2'.'5 for the Ha line (this 
paper) . 



brighter. Thus, despite these crude estimates, we conclude, 
that the absolute brightness of the LBV in PHL 293B is 
one of the highest among the other LBVs and compares 
well with the brig htness of 77 Car, the most lum i nous LBV 
in the Milky Way dMassev 1 . 1200.4 ISmith et al. 1 . 120 lib . All 
this implies that the LBV star in PHL 293B is experiencing 
a strong outburst. 



4.2. Velocities of the stellar wind 

One of the most important features of the VLT/X-shooter 
spectrum is that the P Cygni profiles of the LBV star are 
seen not only in the Balmer hydrogen lines, but also in the 
Paschcn hydrogen and He 1 lines. In particular, no He 1 lines 
with P Cygni profiles were detected and reported in the 
SDSS and VLT/UVES spectra of PHL 293B implying some 
changes in the LBV spectrum after eight years. Examples of 
the lines with the blue absorption minima are shown in Fig. 
H Earlier llzotov fc Thuanl (|2009f > reported the presence of 
He I lines with P Cygni profiles in the spectrum of another 
LBV in the extremely metal-deficient BCD DDO 68. 

In Table [5l we show the terminal velocities Vt of the 
stellar winds for the Balmer, Paschen and He 1 lines and 
the FWHMs of the broad Balmer and Paschen lines of the 
PHL 293B LBV, as derived from the X -shooter (this pa- 
per) and UVES (jlzotov k Thuanl . 120091) spectra. The stel- 
lar wind terminal velocity vt is derived from the wavelength 
difference between the blue absorption minimum and the 
narrow nebular emission line maximum. The v^s derived 
from different lines for two different epochs are similar. 
However, the much larger number of the lines of different 
species in the VLT/X-shooter spectrum allowed us to un- 
cover some new interesting properties. The examination of 
Table [5] shows that the terminal velocities derived from the 
Balmer lines are systematically higher than those derived 
from the Paschen lines and He I lines. Thus, the averaged 
terminal velocities for Balmer, Paschen and He 1 lines are 
818±20, 782±13 and 762±34 km s _1 , respectively, while 
the weighted mean terminal velocities for the same lines 
are 820±8, 783±22 and 764±15 km s" 1 , respectively. These 
differences may have important implications for our under- 
standing of the LBV wind properties. The Balmer absorp- 
tion lines are formed in transitions from the level with the 
principal quantum number n = 2, while for the Paschen 
line transition the lower level is n = 3, which is less popu- 
lated. Therefore, the Balmer absorption lines are formed in 
more outer layers of the expanding envelope as compared 
to the Paschen absorption lines. Probably, the level n = 
2 for He 1 is also less populated as compared to that for 



Table 5. Terminal velocities and FWHMs of the broad 
emission lines in PHL 293B at different epochs 



Line 


X-shooter a 




UVES" 


vt c 


FWHM tr 




FWHM t / 


3750 H12 


833±26 








3770 Hll 


828±26 








3797 H10 


806±24 








3835 H9 


801±27 








3888 H8 


816±24 








3968 H7 


848±24 








4101 RS 


792±24 


878±40 


735 


474±30 


4340 H 7 


836±24 


1006±20 


871 


681±14 


4471 He 1 


787±43 








4861 H/3 


833±25 


966±16 


819 


649±12 


4921 He 1 


697±39 








5016 He 1 


802±30 








5876 He 1 


777±55 








6563 Ha 


787±33 


914±10 


859 


719±10 


6678 He 1 


757±27 








7065 He 1 


752±43 








8753 Pal2 


765±42 








8865 Pall 


799±37 








9017 PalO 


788±54 








9232 Pa9 


775±53 








10941 Pa 7 




1237±160 






12821 Pa/3 




1093±55 







"This 



papc 



llzotov fc Thuanl |2009l ). 

c In km s _1 . 



hydrogen because of its higher energy and lower He abun- 
dance. Therefore, the He 1 absorption lines are formed in 
more inner layers as compared to the Balmer absorption 
lines. Then the differences in the terminal velocities for dif- 
ferent transitions can be explained by a wind accelerating 
outward. 

Another important feature of the VLT/X-shooter 
spectrum is that our derived FWHMs of ~ 1000 km 
s _1 are by a factor of ~ 1.5 higher than those ob- 
tained by llzotov k Thuan (I2009D. but are co nsistent with 



FWHMs derived by llzotov fc Thuanl (|2009l) for another 



low-metallicity LBV in the BCD DDO 68. We also mea- 
sured Vt 's and FWHMs of broad H/3 and Ha emission lines 
in the lower-resolution SDSS spectrum obtained ~ 1 year 
earlier than the VLT/UVES spectrum. We find that a vt of 
~ 800 km s _1 in the SDSS spectrum is consistent with the 
values in Table [5] On the other hand, the FWHMs of the 
broad H/3 and Ha emission lines in the SDSS spectrum, ~ 
1000 km s _1 , arc higher than those in the VLT/UVES spec- 
trum but are similar to the values derived from the VLT/X- 
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shooter spectrum. Despite these differences, all FWHMs 
of the broad lines in the spectra of PHL 293B are signifi- 
cantly broader than those observed in the spectra of high- 
metallicity LBVs whi ch are about 100-200 km s _1 (e.g. 
iLeitherer et al. I . Il994l ). However, they are similar to the 
FWH Ms in some extragala ctic LBVs during their outburst 
phase (|Smith et al. 1 . 120111) . It is possible therefore that the 
LBV in PHL 293B also experiences a strong outburst. 



5. Conclusions 

We have studied here the properties of the narrow-line spec- 
trum of the blue compact dwarf (BCD) galaxy PHL 293B 
and the broad line emission of the luminous blue variable 
(LBV) star in this galaxy from the VLT/X-shooter spec- 
troscopic observations in the wavelength range AA3000 - 
22650A. This dat a has been compared w ith the data ob- 
tained earlier by llzotov fe Thuanl (|2009ft from SDSS and 



VLT/UVES observations. We have arrived at the following 
conclusions: 

1. The oxygen abundance in PHL 293B is 12+log O/H 
= 7.71 ± 0.02, in agr eement with the 12+lo g O/H = 
7.72 ± 0.01 derived by llzotov fc Thuanl (|2009D from the 
VLT/UVE S spectrum an d the 1 2+log O/H = 7.66 ± 0.04 
derived by llzotov et al. I (|2007|) from a lower resolution 
SDSS spectrum. This BCD is the second lowest-metallicity 
galaxy with a detected LBV star, after DDO 68. 

2. We find that the extinction-corrected fluxes of nar- 
row hydrogen lines in the entire wavelength range AA 3000 - 
22650A are in good agreement with the theoretical recom- 
bination values if the extinction coefficient C(H/3) = 0.225, 
derived from the narrow-line Balmer decrement, is adopted 
(Table [I| . This implies that no additional star formation is 
found that is seen in the NIR range but is hidden in the 
visible range. 

3. We detect three v = 1-0 emission lines of 
vibrationally-excitcd molecular hydrogen. The relative 
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fluxes of these lines and non-detection of v = 2-1 and 3- 
1 emission suggest that H2 is collisionally excited in PHL 
293B. 

4. Broad hydrogen emission lines are seen in both the 
optical and NIR spectra of PHL 293B. The broad Ha lumi- 
nosity in the LBV is by a factor of ~ 10 high er than the one 
in the higher-metallicity LBV in NGC 2363 ([Drissen et al. I . 
119971 120011; iPetit et al. I . I2006D and by a factor of ~ 100 
highe r than in the LBV star in DDO 68 ([Izotov k Thuanl . 
[200l . We find some evidence for flux variations of the 
broad lines on a time scale of ~ 8 years. 

5. For the first time we find P Cygni profiles not only 
from Balmer hydrogen lines, but also from Paschen hydro- 
gen and He 1 lines in this LBV. We derive a high terminal 
velocity vt of the stellar wind. The averaged terminal veloc- 
ity from Balmer lines (818±20 km s _1 ) is larger than those 
derived from the Paschen lines (782±13 km s _1 ) and He I 
lines (762±34 km s" 1 ), probably indicating that the wind 
is accelerating outward. These VtS are several times larger 
than the Vt of ~ 100 - 200 km s _1 in some high-metallicity 
counterparts of the Milky Way. However, the high vt is sim- 
ilar to those observe d in some ext r agala ctic LBVs during 
their outburst phase (ISmith et al. I . l20TTt ). 
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